

put /"This is optimal if ifopt = 1 and baseline if ifopt = 0";
put /"ifopt =" ifopt        ;
put "ifdamgis =" ifdamgis   ;
put "ifdamnorm =" ifstandam ;
put /

put "Loopdamcoefgis";
put "loopmeltcoef";
put "loopdisc";
put "loopelast"
put "loopvol";
put "gisvol-1";
put "scc-2015";
put "scc-2050";
put "temp-2050";
put "temp-2100";
put "temp";
put "co2 conc-2050";
put "co2 conc-2100";
put "co2 conc";
put "output-2050";
put "output-2100";
put "emissions-2050";
 put "emissions-2100";
put "emissions-2150";
put "damage fraction-2050";
put "damage fraction-2100";
put "damage fraction-2150";

put "interest rate-2015";
put "interest rate-2100";

put "objective fn";

put "cumulative carbon-2100";
put "cumulative carbon-2200";
put "cumulative carbon-2300";
put "miu-2050";
put "miu-2100";
put "miu-2150";
*put "GISvol-1";
put "GISvol-100";
put "GISvol-200";
put "GISvol-300";
put "GISvol-400";
put "GISvol-500";
put "GISvol-600";
put "GISvol-700";
put "GISvol-800";

put /
loop(kkvol,loop(kkmelt,loop(kkdisc,loop(kkdamgis,

put loopdamcoefgis(kkdamgis);
put loopmeltcoef(kkmelt);
put loopdisc(kkdisc);
put loopelast(kkdisc);
put loopvol(kkvol);
put loopvolgis(kkdamgis,kkmelt,kkdisc,kkvol,'1');

put loopscc(kkdamgis,kkmelt,kkdisc,kkvol,'1');
put loopscc(kkdamgis,kkmelt,kkdisc,kkvol,'2');
put looptemp(kkdamgis,kkmelt,kkdisc,kkvol,'8');
put looptemp(kkdamgis,kkmelt,kkdisc,kkvol,'18');
put looptemp(kkdamgis,kkmelt,kkdisc,kkvol,'38');
put loopmat(kkdamgis,kkmelt,kkdisc,kkvol,'8');
put loopmat(kkdamgis,kkmelt,kkdisc,kkvol,'18');
put loopmat(kkdamgis,kkmelt,kkdisc,kkvol,'38');
put loopy(kkdamgis,kkmelt,kkdisc,kkvol,'8');
put loopy(kkdamgis,kkmelt,kkdisc,kkvol,'18');

put loopem(kkdamgis,kkmelt,kkdisc,kkvol,'8');
put loopem(kkdamgis,kkmelt,kkdisc,kkvol,'18');
put loopem(kkdamgis,kkmelt,kkdisc,kkvol,'28');
put loopdamf(kkdamgis,kkmelt,kkdisc,kkvol,'8');
put loopdamf(kkdamgis,kkmelt,kkdisc,kkvol,'18');
put loopdamf(kkdamgis,kkmelt,kkdisc,kkvol,'28');

put loopr(kkdamgis,kkmelt,kkdisc,kkvol,'1');
put loopr(kkdamgis,kkmelt,kkdisc,kkvol,'18');
put loopobjective(kkdamgis,kkmelt,kkdisc,kkvol);
put loopcca(kkdamgis,kkmelt,kkdisc,kkvol,'18');
put loopcca(kkdamgis,kkmelt,kkdisc,kkvol,'38');
put loopcca(kkdamgis,kkmelt,kkdisc,kkvol,'58');
put loopmiu(kkdamgis,kkmelt,kkdisc,kkvol,'8');
put loopmiu(kkdamgis,kkmelt,kkdisc,kkvol,'18');
put loopmiu(kkdamgis,kkmelt,kkdisc,kkvol,'28');



put loopvolgis(kkdamgis,kkmelt,kkdisc,kkvol,'100');
put loopvolgis(kkdamgis,kkmelt,kkdisc,kkvol,'200');
put loopvolgis(kkdamgis,kkmelt,kkdisc,kkvol,'300');
*put loopvolgis(kkdamgis,kkmelt,kkdisc,kkvol,'400');
*put loopvolgis(kkdamgis,kkmelt,kkdisc,kkvol,'500');
*put loopvolgis(kkdamgis,kkmelt,kkdisc,kkvol,'600');
*put loopvolgis(kkdamgis,kkmelt,kkdisc,kkvol,'700');
*put loopvolgis(kkdamgis,kkmelt,kkdisc,kkvol,'800');
put /
););););

loop(kkvol,loop(kkmelt,loop(kkdisc,loop(kkdamgis,
put /"This is optimal if ifopt = 1 and baseline if ifopt = 0";
put /"ifopt =" ifopt;
put "meltrate =/ " avoldot;
put / "Period" ;
Loop (T, put T.val);
put / "Year" ;
Loop (T, put (2010+(5*T.val) ));

*loop(kkmelt, loop(kkvol, loop(kkdisc, loop(kkdamgis,
put /
put "GISDAMcoef" ;
put loopdamcoefgis(kkdamgis);
put "meltmult" ;
put loopmeltcoef(kkmelt);
put "discrate" ;
put loopdisc(kkdisc);
put "mu elast";
put loopelast(kkdisc);
put "vol lim";
put loopvol(kkvol);

put / "CARBON POLICY" ;
put / "Emissions Control Rate" ;
Loop (T, put loopmiu(kkdamgis,kkmelt,kkdisc,kkvol,t));

put / "Social cost of carbon" ;
Loop (T, put loopscc(kkdamgis,kkmelt,kkdisc,kkvol,t));

put / "Carbon Price (per t CO2)" ;
*Loop (T, put loopmiu(kkdamgis,kkmelt,kkdisc,kkvol,t));
*Loop (T, put cprice.l(T));


put / "ECONOMICS" ;
put / "Output Net Net) " ;
Loop (T, put loopy(kkdamgis,kkmelt,kkdisc,kkvol,t));


*put / "Interest Rate " ;
*Loop (T, put loopr(kkdamgis,kkmelt,kkdisc,kkvol,t));


*put / "Period utilty" ;
*Loop (T, put loopperu(kkdamgis,kkmelt,kkdisc,kkvol,t));

put / "Objective" ;

Loop (T, put loopobjective(kkdamgis,kkmelt,kkdisc,kkvol));
put / "DAMAGE" ;
put / "damages" ;
*Loop (T, put loopmiu(kkdamgis,kkmelt,kkdisc,kkvol,t));
put / "Climate Damages fraction output" ;
Loop (T, put loopdamf(kkdamgis,kkmelt,kkdisc,kkvol,t));
put / "EMISSIONS " ;
put / "Industrial Emissions GTCO2 per year" ;
Loop (T, put loopem(kkdamgis,kkmelt,kkdisc,kkvol,t));

put / "Total Emissions GTCO2 per year" ;
*Loop (T, put E.l(T));
 put / "abatement" ;
*Loop (T, put abatecost.l(t));
 put / "sigma" ;
*Loop (T, put sigma(t));
put / "Cumulative ind emissions" ;
*Loop (T, put cca.l(t));
put / "Cumulative total emissions" ;
Loop (T, put loopcca(kkdamgis,kkmelt,kkdisc,kkvol,t));

put / "CLIMATE" ;
put / "Atmospheric concentrations Gt" ;
*Loop (T, put mat.l(t));
put / "Atmospheric concentrations ppm" ;
Loop (T, put loopmat(kkdamgis,kkmelt,kkdisc,kkvol,t));
put / "Atmospheric Temperature " ;
Loop (T, put looptemp(kkdamgis,kkmelt,kkdisc,kkvol,t));
put / "Forcings" ;
*Loop (T, put forc.l(t));

put / "GREENLAND ICE SHEET" ;
put / "GIS volume, % of original" ;
Loop (T, put loopvolgis(kkdamgis,kkmelt,kkdisc,kkvol,t));
put / "voldot" ;
Loop (T, put loopvdot(kkdamgis,kkmelt,kkdisc,kkvol,t));
put / "tstar*" ;
Loop (T, put looptstar(kkdamgis,kkmelt,kkdisc,kkvol,t));
put / "sign" ;
Loop (T, put loopsign(kkdamgis,kkmelt,kkdisc,kkvol,t));
put / "Global Atmospheric Temperature " ;
Loop (T, put looptemp(kkdamgis,kkmelt,kkdisc,kkvol,t));
put / "Temp diff " ;
Loop (T, put looptd(kkdamgis,kkmelt,kkdisc,kkvol,t));

put /

););););
